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1.  INTRODUCTION 


The  one-dimensional  Navy  Standard  Surf  Model  (NSSM)  has  been  shown  to  be 
very  robust,  but  it  has  its  limitations.  Since  NSSM  assumes  parallel  bottom  contours  in 
the  surf  zone  (generally  depths  of  8  m  or  less),  it  cannot  account  for  longshore  variations 
of  bathymetry  or  forcing.  NSSM  can  produce  inaccurate  wave  and  longshore  current 
estimations  for  areas  with  complicated  bathymetry  where  a  2D  nearshore  hydrodynamic 
model  is  needed.  The  DelftSD  modeling  system,  developed  by  Delft  Hydraulics 
(www.wldelft.nl),  is  capable  of  simulating  hydrodynamic  processes  due  to  waves,  tides, 
rivers,  winds  and  coastal  currents.  Using  a  version  of  DelftSD  with  roller  dynamics, 
Morris  (2001)  has  shown  that  DelftSD  produces  good  results  when  compared  with 
DELILAH  and  Duck94  data  at  Duck,  North  Carolina,  and  Torrey  Pines  and  Santa 
Barbara  data  in  California.  The  major  remaining  difficulty  is  the  specification  of  bottom 
friction.  Morris  reported  that  different  values  of  roughness  for  bottom  friction  should  be 
used  depending  on  whether  the  beach  in  the  domain  is  barred  or  planar.  But,  there  is  no 
established  guideline  to  choose  these  values.  Before  making  DelftSD  an  operational 
model  for  surf  applications,  this  remaining  bottom  friction  issue  needs  to  be  further 
examined. 

Improvements  of  the  wave/roller  model  and  radiation  at  the  side  boundaries  have 
been  implemented  in  DelftSD  in  recent  years  (Roelvink,  2003).  The  main  objective  of 
this  report  is  to  evaluate  the  performance  of  the  new  DelftSD  and  its  sensitivity  to  bottom 
friction  selection.  In  our  study,  the  October  data  of  Duck94  is  used  for  evaluation  and 
calibration.  The  Santa  Barbra  data  is  used  for  validation,  verifying  the  various  model 
parameter  selections  based  on  Duck94  data. 

2.  MODEL  DESCRIPTIONS 

The  DelftSD  system  uses  two  main  modules  for  simulating  nearshore  wave- 
induced  hydrodynamic  processes.  The  WAVE  module  includes  the  latest  version  of 
SWAN  model  for  computing  propagation  and  generation  of  waves.  Hydrodynamic  flow 
is  simulated  with  the  LLOW  module  (WL  Delft  Hydraulics,  2001),  which  solves  the 
unsteady  shallow  water  equations  in  two  (depth- averaged)  or  three  dimensions.  Lor 
applications  with  dominant  wave-induced  flow,  only  the  two-dimensional  mode  is  used. 
The  model  can  be  mn  with  both  one-way  forcing  or  with  feedback  between  the  two 
modules. 

DelftSD  can  be  run  in  Cartesian  (equidistant  or  stretched)  or  curvilinear 
coordinates;  all  necessary  grid  generation  software  for  creating  curvilinear  grids  is 
included  with  the  DelftSD  package.  Although  a  highly  flexible  tool  for  various 
applications,  this  component  of  the  nearshore  modeling  system  was  best  suited  for  a 
domain  that  would  extend  from  the  shoreline  to  about  one  kilometer  seaward  (Dykes  et 
ah,  2003). 

Previous  field  data  shows  that  the  peak  longshore  current  occurred  in  the  trough 
whereas  the  model  peak  occurred  at  the  beginning  of  the  bar  (Visser,  1982).  The  roller 
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formulation  (e.g.,  Ruessink,  et  al.,  2001)  was  proposed  to  delay  of  momentum  release 
from  the  wave  breaking.  So  the  inclusion  of  this  formulation  was  mainly  designed  to 
move  the  velocity  peak  towards  the  trough  of  the  bar.  In  the  DelftSD  roller 
implementation  manual,  Roelvink  (2003)  described  the  roller  model  as  follows: 

•  It  solves  a  wave  energy  balance  and  roller  energy  balance  within  the  FLOW  model,  using 
a  trimmed-down  version  of  the  ‘DIFU’  advection-diffusion  solver. 

•  Based  on  the  wave  energy  and  roller  energy,  radiation  stresses  and  gradients  of  these 
stresses  are  computed  that  replace  the  conventional  wave  forces  as  derived  from  the 
WAVE  model 

•  The  wave  direction  is  still  derived  from  an  initial  WAVE  run. 

The  full  equations  for  the  flow  and  roller  model  are  not  repeated  here  (Roelvink,  2003; 
Reniers  et  ah,  2004). 

3.  FIELD  DATA  AND  MODEL  SETUP 
3.1  Duck94  Cases 

Duck94  field  experiment  was  conducted  in  August  and  October  near  the  Army  Corps 
of  Engineers’  Field  Research  Facility  pier  located  in  Duck,  North  Carolina.  The 
instmment  layout  is  shown  in  Fig.  1.  The  so  called  mini-grid  is  marked  by  the  box.  A 
bathymetry  survey  for  the  mini-grid  was  conducted  daily  during  the  intensive  study 
period  in  October.  The  wave,  wind  and  tide  conditions  during  October  were  shown  in 
Fig.  2.  A  storm  from  the  southeast  occurred  around  October  15  with  a  peak  significant 
wave  height  reaching  4  m.  The  beach  profile  was  significantly  changed  after  the  storm  as 
illustrated  in  Fig.  3. 


(  from  www.frf.usace.armv.mil/duck94/duck94  overview. stm) 
Fig.  1-  Instmment  layout  at  Duck94. 
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( from  www.frf.usace.armv.mil/duck94/duck94  overview.stm) 
Fig.  2  -  Environmental  conditions  during  October 


Figure  5.  Profile  change  aiong  the  Elgar  &  Guza  cross-shore  array  during  DUCK94 


( from  www.frf.usace.armv.mil/duck94/duck94  overview.stm) 
Fig.  3  -  Profile  changes  before  and  after  the  October  15  storm. 
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The  WAVE  model  grid  consisted  of  86  by  131  grid  points  with  grid  spacing  of 
10  m  in  the  cross-shore  direction  and  15  m  along  the  shore.  The  offshore  boundary  was 
located  at  the  8  m  directional  wave  gage  array.  The  smaller  FLOW  grid  lay  inside  the 
WAVE  domain  with  86  by  117  grid  points  with  grid  spacing  of  10  m  in  the  cross-shore 
direction  and  15  m  along  the  shore.  The  dates  (October  5-21  and  26)  chosen  for  model 
validation  was  based  on  the  availability  of  daily  mini-grid  bathymetry  surveys.  A 
bathymetry  survey  for  a  much  larger  area  at  a  coarser  spacing  was  conducted  every  other 
month.  By  examining  the  coarser  September  and  December  data,  there  is  little 
bathymetry  change  beyond  the  5-m  contour.  Most  of  the  depth  changes  are  within  the 
mini-grid  region.  To  build  the  model  input  depth  files,  the  coarser  data  was  merged  with 
daily  mini-grid  data.  Because  of  the  October  15  storm,  September  data  was  used  for 
merging  for  those  dates  before  October  14  and  November  data  was  used  afterwards. 

3.2  Santa  Barbara  Cases 

There  are  four  cases  from  NSTS  field  experiments  at  Leadbetter  Beach,  Santa 
Barbara,  California  in  Febmary,  1980.  The  offshore  waves  consisted  mainly  of  swells 
with  narrow  banded  frequency  and  direction.  The  beach  was  planar  with  a  steeper  slope 
than  at  Duck.  The  nearshore  depth  contours  during  the  experiment  were  nearly  straight 
and  parallel  to  the  beach.  The  2D  bathymetry  for  model  input  was  constmcted  based  on  a 
single  trace  of  the  beach  profile.  The  wave  grid  was  54  by  71  grids  with  2  m  in  x  and  5 
m  in  y  grid  size.  The  flow  grid  is  smaller  than  wave  grid  with  52  by  41  at  the  same  grid 
resolution. 

4.  MODEL  RESULTS 
4.1  Duck94  Cases 

Although  model  mns  are  made  from  October  5  through  21  and  26  at  every  3 
hours,  only  sample  results  are  presented  for  discussion.  Both  low-  and  high-tides  cases  as 
well  as  cases  when  waves  came  from  both  northeast  and  southeast  were  chosen.  Cases 
were  also  chosen  based  on  the  occurrence  of  storms,  i.e.  before  and  after  the  October  15 
storm. 


4.1.1  Roller  vs.  No  Roller 

As  mentioned  in  the  introduction,  the  purpose  of  roller  implementation  is  to  move 
the  velocity  peak  further  towards  the  trough  of  the  bar.  Comparisons  between  cases  with 
and  without  roller  are  shown  in  Fig  4.  The  top  panel  shows  the  beach  profile  and  gage 
locations.  The  middle  panel  shows  the  significant  wave  height  comparison  and  third 
panel  shows  the  longshore  current  comparison.  With  the  roller  turned  on,  wave  height  is 
computed  by  the  Delft3D  roller/wave  module  and  only  uses  SWAN  computed  wave 
angles.  The  plotted  wave  height  is  derived  from  the  wave  energy  term  in  the  binary 
“trim”  (map)  file.  With  roller  off,  the  plotted  wave  height  computed  from  SWAN  is 
derived  from  the  binary  “com”  (communication)  file.  In  all  cases,  the  inclusion  of  roller 
improves  the  agreement  by  moving  the  velocity  peak  towards  shore  as  expected. 
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(d)  Oct.  26,  EST  2200 

Eig.  4  -  Comparisons  between  cases  with  and  without  roller.  EST  represents  Eastern 
Standard  Time. 

4.1.2  Breaking  Dissipation  Formulation 

In  DelftSD,  the  use  of  new  features  which  have  not  been  made  to  be  a  standard 
version  is  specified  under  the  data  group  “additional  parameters”.  Keyword  and  its 
corresponding  value  are  entered.  One  available  keyword  is  Gamdis  for  specifying  the 
gamma  value,  i.e.  the  significant  wave  height  to  depth  ratio  associated  with  depth- 
induced  breaking.  In  the  present  version,  if  one  uses  the  stationary  roller  model,  the 
Baldock  dissipation  model  is  used,  for  which  Gamdis  is  set  at  the  default  value  of  0.55. 
If  one  sets  the  Gamdis  keyword  at  -1,  variable  gamma  value  is  used  using  Ruessink's 
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formulation  (Ruessink  et  al.,  2003).  The  comparisons  between  Baldock  and  Ruessink 
breaking  dissipation  formulation  are  shown  in  Fig.  5.  In  general,  the  Ruessink’ s 
formulation  produces  better  agreement  in  longshore  current  with  a  narrower  velocity 
distribution  profile,  i.e.  smaller  offshore  values  and  more  peaked  at  velocity  maximum. 
The  shape  of  the  velocity  profile  also  depends  on  the  empirical  constant  of  turbulent  eddy 
viscosity  associated  with  mixing  (Battjes,  1975).  The  default  value  is  1.  The  same  value 
is  also  used  by  NSSM  (Hsu,  et  al,  2000).  Note  that  the  background  horizontal  eddy 
viscosity  as  defined  in  the  data  group  “physical  parameters”  is  not  used  when  roller  is 
turned  on. 
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(a)  Oct.  10,  EST  1600 
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(c)  Oct.  19,  EST  1000 


DUCK94-9410262200 


(d)  Oct.  26,  EST  2200 

Eig.  5  -  Comparison  of  breaking  dissipation  formulations. 

4.1.3  Breaker  Delay 

An  optional  feature  called  'breaker  delay’  is  available  through  the  keyword  E_lam 
with  a  value  -2,  the  description  is  given  in  Roelvink  et  al  (2000).  It  applies  a  weighted 
averaging  of  the  depth  used  in  the  wave  dissipation  model,  to  a  number  of  wavelengths 
seaward.  The  comparisons  between  cases  with  and  without  breaker  delay  are  shown  in 
Eig.  6.  Except  for  Oct.  26,  all  other  cases  show  a  slight  movement  of  peak  velocity 
towards  trough.  Eor  case  of  Oct.  26,  the  improvement  is  negligible. 
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DUCK94-9410191000 


(c)  Oct.  19,  EST  1000 


DUCK94-9410262200 


(d)  Oct.  26,  EST  2200 

Eig.  6  -  Comparison  between  cases  with  and  without  applying  breaker  delay 

4.1.4  Turn  off  the  Roller  Stress  at  very  shallow  water 

In  very  shallow  water,  say  depth  less  than  0.5  m,  DelftSD  occasionally  generates 
large  currents,  due  to  unrealistic  roller  or  wave  forces.  In  DelftSD,  roller  stress  is  turned 
off  for  depth  shallower  than  2  times  the  speeified  threshold  depth,  which  is  used  to  define 
the  depth  above  which  a  grid  cell  is  considered  to  be  wet.  Designated  as  dryflc  under  the 
data  group  “numerical  parameters”,  its  default  value  is  0. 1  m.  Erom  our  experience,  the 
value  of  dryflc  needs  to  be  slightly  increased  to  avoid  unrealistieally  large  eurrent  at 
shallow  depth.  The  impaet  of  the  seleetion  of  dryflc  is  illustrated  in  Eig.  7.  The 
significant  reduction  of  the  second  peak  near  shoreline  at  dryflc  =  0.2  m,  i.e.  roller  stress 
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turned  off  at  0.4  m,  is  very  clear  in  case  (a).  In  the  other  cases,  the  reduction  is  small.  It 
is  important  to  note  that  the  selection  of  dryflc  does  not  have  much  impact  to  the  rest  of 
the  velocity  distribution  away  from  the  more  shallow  water. 


(a)  Oct.  10,  EST  1600 


(b)  Oct.  12,  EST  1300 


11 


(c)  Oct.  19,  EST  1000 


(d)  Oct.  26,  EST  2200 

Eig.  7  -  The  effects  of  turning  off  roller  stress  at  very  shallow  water. 

4.1.5  The  Effect  of  Bottom  Friction  Selection 


Eor  2-D  simulations,  three  roughness  formulations,  i.e.  Manning,  White- 
Colebrook  (W-C)  and  Chezy,  can  be  selected  at  the  roughness  menu  under  the  data  group 
“physical  parameters”.  By  definition,  the  bottom  friction  coefficient,  C/,  is  related  to 
Chezy  roughness  coefficient  C  by 

Cf=glC^  (1) 

where  g  is  the  acceleration  due  to  gravity.  The  default  Chezy  formulation  is  for 
specifying  a  constant  Chezy  value  for  the  whole  domain,  i.e.  no  depth  dependency.  Eor 
Manning  formulation,  the  Chezy  coefficient  is  related  to  the  Manning  roughness 
coefficient,  n,  as. 
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C  =  h''^  In 


(2) 


where  h  is  the  water  depth.  For  W-C  formulation,  one  needs  to  specify  the  equivalent 
geometrical  roughness  of  Nikuradse,  ks.  The  Chezy  coefficient  is  calculated  from: 

C  =  18-log(l2-/i/kJ  (3) 

In  Morris’s  investigation  (2001),  he  chose  the  W-C  formulation.  The  optimal 
value  of  ks  based  on  Duck  runs  was  found  to  be  0.003  m  for  the  barred  beach 
corresponding  cross-shore  averaged  C/  about  0.002.  He  also  reported  that  the  optimal 
value  of  ks  based  on  planar  beach  at  Torrey  Pines  and  Santa  Barbra  beaches  was  0.009  m 
with  a  corresponding  averaged  C/  at  0.004.  As  for  the  Manning  formulation,  the 
suggested  default  value  for  n  is  0.02  in  the  Delft3D  FLOW  manual.  Gmnnet  et  al.  (2004) 
uses  a  value  of  0.0225  for  a  barrier  island  sediment  study  using  Delft3D.  The  default 
value  for  Chezy  formulation  is  65,  equivalent  to  a  bottom  friction  coefficient  of  0.0023. 
The  depth  dependence  of  both  Manning  and  W-C  formulations  is  illustrated  in  Fig.  8  for 
a  linear  beach.  Both  formulations  show  a  weak  dependence  of  depth  until  depth  is  very 
shallow,  say  0.5  m.  It  is  noted  that  the  curve  with  Manning  n  at  0.02  and  W-C  ks  at  0.009 
are  very  similar.  The  averaged  C/  value  will  depend  on  the  shape  and  depth  of  the  real 
beach  profile.  For  the  linear  beach  in  Fig.  8,  the  corresponding  cross-shore  averaged  C/ 
value  for  n  =  2  is  about  0.0028,  and  is  about  0.0024  for  ks  at  0.009. 


Fig.  8  -  Depth  dependence  of  bottom  friction  coefficient,  Cf ,  for  Manning  and  White- 
Colebrook  roughness  formulations. 
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To  evaluate  the  effects  of  bottom  friction  dependence,  comparison  plots  are 
presented  in  Fig.  9.  The  figures  show  that  Morris’  selection  of  W-C  at  ks  =  0.003  based 
on  DELILAH  field  data  still  works  reasonably  well  for  Duck94.  The  cases  using  default 
Chezy  C  at  65  tend  to  over-predict  the  velocity  peaks  whereas  the  cases  using  C  at  55 
(equivalent  to  C/  at  0.0032)  match  the  data  slightly  better.  The  Manning  formulation  at 
n  =  0.02  works  well. 


DUCK94  -9410101600 


(a)  Oct.  10,  EST  1600 


DUCK94  -9410121300 


(b)  Oct.  12,  EST  1300 
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DUCKP4--9410191000 


(c)  Oct.  19,  EST  1000 


DUCK94--9410262200 


Crossshore  Distance  (meters) 

(d)  Oct.  26,  EST  2200 

Eig.  9  -  The  comparison  of  bottom  friction  selections 

To  study  the  impact  of  bottom  friction  further,  64  out  of  the  computed  116  cases 
of  Duck94  data  were  used  for  statistical  computation.  The  cases  with  rip  current  or  big 
eddies  are  not  appropriate  for  evaluating  the  effect  of  bottom  friction,  and  therefore  were 
not  included.  The  eddy  formation  is  associated  with  the  instability  at  the  southern 
boundary  when  waves  are  coming  from  southeast.  Eor  Neumann  boundary  conditions  to 
work  well,  the  bathymetry  longshore  near  the  side  boundaries  of  SWAN  needs  to  be 
uniform.  More  than  20  grid  spaces  of  alongshore  uniform  depth  values  are  provided  for 
the  northern  boundary  whereas  only  5  were  provided  in  the  southern  boundary.  It  turns 
out  that  5  grid  spaces  are  not  enough  for  Neumann  boundary  condition  to  work  well.  As 
a  rule  of  thumb,  more  than  15  grid  spaces  should  be  used  in  both  boundaries  in  any  future 
setup.  It  is  useful  to  set  a  few  side  boundary  grid  cells  in  the  flow  bathymetry  to  the  same 
depth  as  SWAN  to  achieve  further  smooth  transition. 
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In  the  statistics  computation,  all  cases  are  run  with  the  following  conditions: 

1.  roller  and  breaker  delay  on, 

2.  breaking  dissipation  formulation  by  Ruessink  et  al.  (2003), 

3.  roller  stress  turned  off  at  0.4  m, 

4.  wind  turned  on, 

5.  Neumann  boundary  conditions  on  side  boundaries,  and 

6.  input  from  the  directional  spectra  from  the  8-m  gage  array. 

Since  wave  height  prediction  has  been  shown  to  be  accurate,  no  more  discussion  is 
presented  here.  Just  the  scatter  plot  is  presented  in  Fig.  10a  for  reference.  With  the 
default  of  bottom  friction  selection  C  at  65,  the  scatter  plot  for  longshore  current  is  shown 
in  Fig.  10b.  For  wave  height,  the  agreement  is  very  good  with  a  correlation  coefficient  at 
0.98  and  RMS  error  at  0.14  m.  For  longshore  current,  it  is  evident  that  friction  setting  at 
C  =  65  is  too  low  as  shown  by  the  slope  of  the  linear  regression  line.  Additional  runs 
were  made  with  C  =  55,  equivalent  to  a  constant  C/=  0.0032.  The  resulting  scatter  plot  is 
shown  in  Fig.  1 1  shows  an  improvement  of  agreement.  Additional  mns  were  made  with 
Ks  =  0.009  and  n  =  0.002  and  are  presented  in  Figs.  12  and  13. 


Wave  Height 


(a)  Wave  Height  Comparison 
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Longshore  Current 


Measured  (m/s) 


(b)  Longshore  Current  Comparison 


Fig.  10  -  Scatter  plots  for  (a)  wave  height  and  (b)  longshore  current  using  default  Chezy 
formulation  with  C  =  65.  R  is  linear  correlation  coefficient;  Slope  is  the  slope  of  the 
linear  regression  line  (solid  line).  N  is  the  number  of  observations. 


Longshore  Current 


Fig.  1 1  -  Scatter  plot  for  longshore  current  excluding  cases  with  eddy  formation  using 
Chezy  with  C  =  55. 
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Longshore  Current 


Measured  (m/s) 


Fig.  12  -  Scatter  plot  for  longshore  current  excluding  cases  with  eddy  formation  using  W- 
C  formulation  at  ks  =  0.003. 


Longshore  Current 


Fig.  13  -  Scatter  plot  for  longshore  current  excluding  cases  with  eddy  formation  using 
Manning  formulation  at  n  =  0.002. 


The  skill  statistics  of  four  different  bottom  friction  formulations  are  summarized  in 
Table  1.  In  addition  to  root  mean  squares  error  (RMSE),  mean  absolute  gross  error 
(MAGE)  and  mean  bias  error  (MBE)  are  also  listed.  The  model  performance  at  default 
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Chezy  C  =  65  (equivalent  to  Cf  at  0.0023)  is  worse  that  C  =55  (equivalent  to  C/  at 
0.0032),  indicating  the  default  bottom  friction  is  slightly  too  low.  This  is  consistent  with 
previous  field  measurements  of  bottom  friction  at  Duck  (Whitford  and  Thornton,  1996). 
Overall,  the  performance  of  all  three  different  formulations  is  very  similar.  The  Manning 
formulation  with  n  =  0.02  and  Chezy  formulation  with  C=55  produce  better  performance. 

Table  1  -  Skill  statistics  for  different  bottom  friction  selections 


RMSE(m/s) 

MAGE  (m/s) 

MBE  (m/s) 

R 

Slope 

N 

Chezy  C=65 

0.18 

0.13 

0.024 

0.914 

1.12 

533 

Chezy  C  =55 

0.14 

0.10 

-0.009 

0.918 

0.9 

533 

W-C  k  =0.003 

0.20 

0.15 

0.041 

0.912 

1.20 

533 

Manning  n=0.02 

0.14 

0.11 

0.002 

0.919 

0.95 

533 

The  statistics  of  RMS  error  is  influenced  by  the  large  number  of  low  velocity  data  points. 
Fig.  14  shows  the  scatter  plot  under  Chezy  at  C  =  55  where  no  measured  data  under  0. 1  m 
are  included.  The  RMS  value  is  increased  as  expected.  The  skill  statistics  of  data 
excluded  velocity  below  0.1  m/s  is  shown  in  Table  2. 


Longshore  current 


Fig.  14  -  Scatter  plot  for  longshore  current  excluding  measured  velocity  less  than  0.1 
m/s  using  Manning  formulation  at  n  =  0.02. 
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Table  2  -  Skill  statistics  for  different  bottom  friction  selections  excluding  measured 
velocity  less  than  0.1  m/s. 


RMSE(m) 

MAGE  (m/s) 

MBE  (m/s) 

R 

Slope 

N 

Chezy  C=65 

0.2 

0.16 

0.052 

0.914 

1.09 

374 

Chezy  C  =55 

0.16 

0.12 

0.001 

0.918 

0.87 

374 

W-C  k,  =0.003 

0.23 

0.18 

0.076 

0.91 

1.17 

374 

Manning  n=0.02 

0.16 

0.13 

0.016 

0.917 

0.93 

374 

4.2  Santa  Barbara 

Unlike  the  barred  beach  at  Duck,  Santa  Barbara  has  planar  and  steeper  beach,  and 
therefore  is  useful  in  testing  the  robustness  of  DelftSD  features.  Since  the  four  available 
cases  are  similar  in  wave  input  and  beach  profile,  only  one  case,  i.e.  Feb.  03,  1980,  is 
presented. 

4.2.1  Roller  vs.  No  Roller 

Comparison  of  results  with  and  without  roller  is  shown  in  Fig.  15.  Similarly  to  Duck, 
the  case  with  roller  gives  better  longshore  current  results.  The  sharp  decay  of  the 
longshore  current  at  the  beach  for  the  roller  case  is  partly  because  the  roller  stress  was 
turned  off  at  0.4  m.  The  case  was  mn  with  the  Manning  formulation  at  n  =  0.015. 


NSTS  Santa  Barbara  Feb.  03, 1980 


Fig.  15  -  Comparison  between  with  and  without  roller. 
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4.2.2  Breaking  Dissipation  formulation 

The  comparison  of  default  and  Ruessink  formulations  is  presented  in  Fig.  16. 
The  Ruessink  approach  also  gives  better  results  for  a  planar  beach. 


NSTS  Santa  Barbara  Feb.  03, 1980 


Fig.  16-  Comparison  between  different  breaking  dissipation  formulations. 

4.2.3  Breaker  Delay 

The  comparison  between  that  which  was  with  and  without  breaker  delay  is 
presented  in  Fig.  17  The  breaker  delay  produces  worse  agreement  in  longshore  current. 


NSTS  Santa  Barbara  Feb.  03, 1980 


Fig.  17-  Comparison  between  with  and  without  breaker  delay. 


21 


4.2.4  The  Effect  of  Bottom  Friction  Selection 


The  comparison  between  four  different  bottom  friction  selections,  Chezy  with  C  = 
65,  C  =  55  and  Manning  with  n  =  0.0.015  and  n  =  0.02,  is  presented  in  Fig.  18  Compared 
with  Duck  results,  it  is  evident  that  sensitivity  to  the  empirical  bottom  friction  value  is 
much  higher  at  the  steeper  Santa  Barbara  beach.  Unlike  the  Duck94  cases,  the  default 
choice  of  Chezy  with  C  =  65  (C/=  0.0023)  produces  better  agreement  than  C  =  55  (C/  = 
0.0032).  In  Morris’  investigation,  the  optimal  value  (ks  =  0.009  corresponding  averaged 
Cf  at  0.004)  required  at  Santa  Barbara  is  twice  higher  than  at  Duck  {ks  =  0.003 
corresponding  averaged  Cf  at  0.002).  So  our  trend  is  reverse  to  Morris  findings.  It  is 
difficult  to  isolate  what  caused  the  difference  between  these  two  different  versions  of 
Delft3D,  because  a  lot  of  improvements  have  been  implemented.  One  reason  for  the 
difference  may  be  associated  with  the  different  breaking  dissipation  formulation  used. 


Fig.  18-  Comparison  of  Delft3D  results  between  different  bottom  friction  selections. 

5.  SUMMARY  AND  CONCLUSIONS 

The  main  focus  of  this  investigation  is  to  examine  the  effects  of  bottom  frictions 
selection.  Both  Duck94  and  NSTS  Santa  Barbara  data  were  used  for  evaluating  Delft3D 
performance.  All  three  bottom  friction  formulations,  i.e.  Chezy,  White-Colebrook  and 
Manning,  were  evaluated  and  all  produced  good  longshore  current  results,  if  proper 
empirical  constant  is  used.  The  RMS  error  for  longshore  current  is  about  0.2  m/s.  Over 
barred  beach  at  Duck,  the  optimal  values  are  n  =  0.02  (Manning)  and  C  =  55  (Chezy,  Cf 
=  0.0032).  Over  the  steeper  planar  beach  at  Santa  Barbara,  the  optimal  values  are 
Manning  with  n  =  0.015  and  Chezy  with  C  =  65  (C/=  0.0023).  The  difference  in  optimal 
values  (i.e.  C/  from  0.0023  to  0.0032)  between  different  beaches  is  therefore  not  as  high 
as  indicated  by  a  previous  study  (i.e.  Cf  from  0.002  to  0.004,  Morris,  2003).  This 
difference  may  be  attributed  to  the  difference  in  breaking  dissipation  formulations  used. 
Further  evaluation/validation  for  some  other  beaches  is  needed  to  provide  a  general 
guideline  for  the  selecting  optimal  values. 

Many  other  Delft3D  model  parameters  and  options  were  also  evaluated.  The 
inclusion  of  roller  is  shown  to  improve  the  longshore  current  prediction  for  both  Duck 
and  Santa  Barbara  beaches.  As  the  selection  of  breaking  dissipation  formulation,  the 
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Ruessink  formulation  with  variable  gamma  value  (significant  wave  height  to  depth  ratio) 
provides  better  results  than  the  default  setting.  The  optional  breaker  delay  slightly 
improves  the  longshore  current  prediction  at  Duck,  but  it  produces  a  worse  result  at  Santa 
Barbara  beach.  To  avoid  spuriously  high  currents  at  very  shallow  depth,  roller  stress  is 
recommended  to  be  turned  off  at  shallow  depth,  say  0.4  m.  In  general,  DelftSD  has 
shown  to  be  robust  and  accurate  in  predicting  nearshore  flows. 
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